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ABSTRACT 

We present a basic model to understand the redshift evolution of the UV luminosity and stellar 
mass functions of Lyman Break Galaxies. Our minimal approach is based on the assumption that 
the luminosity and stellar mass of a galaxy is related to its dark matter halo assembly and gas infall 
rate. Specifically, galaxies experience a burst of star formation at the halo assembly time, followed 
by a constant star formation rate, representing a secular star formation activity sustained by steady 
gas accretion. The model is calibrated by constructing a galaxy luminosity versus halo mass relation 
at z = 4 via abudance matching. After this luminosity calibration, our model naturally fits the 
z = 4 stellar mass function. Furthermore, it correctly predicts the evolution of both luminosity and 
stellar mass functions over 13 billion years, from z — to z — 8. Star formation from steady gas 
accretion is the dominant contribution to the galaxy luminosity at all redshifts, but the initial burst 
of star formation contributes significantly to the stellar mass assembly and evolution of the stellar 
mass function. While the details of star formation efficiency and feedback are hidden within our 
calibrated luminosity versus halo mass relation, our study highlights that the primary driver of galaxy 
evolution across cosmic time is the build-up of dark matter halos, without the need to invoke a redshift 
dependent efficiency in converting gas into stars. 

Subject headings: cosmology: theory — galaxies: high-redshift — stars: formation 



1. INTRODUCTION 

The galaxy luminosity function (LF) and the stel- 
lar mass function (MF), along with their redshift evo- 
lution, summarize key information on galaxy proper- 
ties and on their evolution with cosmic time. The 
rest-frame UV 1500 A LF in particular can be rela- 
tively easily traced with current technol ogy across the 
whol e redshift range from z ~ (e.g. , iBlanton et al.l 
120031: lArnouts etall 120051: lOesch et all 120101) to z~^ 
10 fe.gTlBouwens et all 120071: IReddv fc SteideH [200l 
Bradl ey et al.l 1201% lOesch et al.l 1201% IZheng et all 
2012f ). the current frontier of detection of Lyman Break 
Galaxy (LBG) populations. Similarly, the stellar MF can 
be derived from obs e rvations in the rest-frame optical 
(lArnouts et all 120071: IStark et all 1201% IGonzalez et all 
120111 ). These data enable a self-consistent comparison 
of star forming galaxies over the entire span of cos- 
mic history. In a A Cold Dark Matter (ACDM) cos- 
mology, a necessary step is to link the properties of 
galaxies to those of their host dark-matter (DM) ha- 
los. This can be d one through numeric a l simulations 
(ISmith et al.l 120111: iFinlator et all 120111: IJaacks et all 
120121 ), semi-analytic models based on DM halo merge r 
trees (|Somerville fc Primackl [l999t iCroton et all |2006), 
as well as halo occupation distribution models. The 
latter give the p robability that a halo of mass My, 
hosts a galaxy (iJing. et al. I 119981 : IPeacock fc Smith! 
2000), which can be generalized into a conditional lu- 
minosity function (CLF) modeling, giving the proba- 
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bility a halo of mass Mh hosts a galaxy with luminos- 
ity L (Ivan den Bosch et"alll20rjl iVale fc Ostrikerl 12004 
ICooravll2005T) ~ 

There are pros and cons with every method. Simula- 
tions enable a self-consistent treatment of the galaxies in 
a cosmological framework, but use several assumptions 
regarding "sub-grid" physics regulating star-formation 
and feedback processes, on scales smaller than the nu- 
merical resolution. Semi-analytic models allow for the 
exploration of a large volume of parameter space, but rely 
on tens of free parameters in order to reproduce observed 
galaxy properties. The CLF approach provides a L(Mh) 
relation between galaxy luminosity and DM halo mass, 
deriv ed at each redshift through "abundance ma tching" 
(e.g.. lMo fc Fukugitalfl996HVale fc Ostrikeil200l , which 
generally includes a duty cycle parameter so as to pop- 
ulate with UV lum i nous galaxies o nly a fraction of DM 
halos (|Cooravll200llLee et al.|[2009h . The CLF approach 
is quite successful in providing a description of the LF, 
but it does not provide a physical explanation for it. 

Our approach is to aim at identifying the key drivers of 
the evolution of galaxy properties with the l east amount 
of ass umptions. We started by showing in iTrenti et all 
(2010) that the LF at z > 5 can be successfully mod- 
eled by assuming that UV bright galaxies are present, 
at any cosmic epoch, only in halos assembled within At 
(At ~ 10 8 Myr). The outcome of this assumption is a 
duty-cycle which is physically motivated, defined with- 
out free p arameters, and depe ndent on redshift and halo 
mass. The lTrenti et all ()2010D model well reproduces the 
observed rest-frame UV LF evolution out to the largest 
redshifts currently probed by observations (z ~ 10). Un- 
derstandably, however, it can not be extrapolated down 
to redshifts z < 4, since it is to be expected that, at such 
late epochs, DM halos older than a few 10 s yr are likely 
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Fig. 1. — Observed UV luminosity functions (LFs) at low-z (left) and high-z (right) with datapoints from different groups (see labels 
and Table[T]l. Observations are quantitatively described by our LF model at all redshifts. In particular, we reproduce the evolution of the 
faint-end slope a, from its shallow low-z value a(z ~ 0) ~ —1.3 to the steepening at z > 6 where a < —1.7. Shaded areas represent 68% 
confidence regions at z = 4 and z = 6. 



to host UV bright galaxies. 

In this Letter we expand the lTrenti et al.l (|2010D model 
by making the more realistic assumption that, at any 
epoch, all massive DM halos host a galaxy with a star 
formation history (SFH) that is related to the time of 
halo assembly. We anchor our model to the observed 
LF at z = 4, and evolve it towards higher (z « 8) and 
lower (z r; 0) redshifts with a simple physical prescrip- 
tion which enables us to explore the origin of the ob- 
served UV LF evolution. Our new model features (i) a 
burst of star formation at halo assembly time, followed 
by (ii) constant star formation with rate inversely pro- 
portional to the halo assembly time (halos at a given 
mass and different redshift accrete the same gas but over 
a different timescale). These assumptions, calibrated at 
2 = 4, are able to reproduce the evolution of the LF, 
cosmic mass density (pm* ) an d specific star formation 
rate (sSFR) across 13 billion years of cosmic time. This 
good match between model and observations is achieved 
with a dominant contribution to the UV luminosity, at 
all epochs, of a continuous mode of star formation, fu- 
eled by gas accretion. However, the initial burst of star 
formation contributes similarly to the stellar mass build 
up of galaxies. 

This Letter is organized as follows. Section [2] sum- 
marizes the observational datasets we aim at modeling. 
Section [3] describes the model and its calibration. Sec- 
tion |4] presents our results and discusses model uncer- 
tainty. We highlight some concluding remarks in Sec- 
tion El Wc adopt WMAP5 cosmology: tt A = 0.72, 
Q m0 = 0.28, n ft .n = 0.0462 , a 8 = 0.817, n s = 0.96, 
h = 0.7 (jKomatsu et al.ll2009fh 



2. OBSERVATIONAL DATA 



Figure [T] and Table [T] summarize the literature com- 
pilation of observed UV LFs we use, including the pa- 
rameters for the best-fit Schechter LF functions (<j>(L) = 
<jf{L/L*) a exp — (L/L*)). Observed stellar mass densi- 
ties (pm.) integrated above M* )m j„ = 10 8 M Q and spe- 
cific star formation rate at M* = 5 x 10 9 M Q (sSFR) are 
shown in Figure O The figure also includes the cosmic 
star forma tion rate dens i ty , obtained by converting, 
using the iMadau et al.1 (|1998h relation, the luminosity 
density p^ integrated to L TO , n = 0.05 L*_ s (correspond- 
ing to Mab — —17.7). Over plotted to observations are 
our model predictions, obtained with the prescriptions 
described below. 

3. MODEL DESCRIPTION 

Our model links the UV LF to abundance of DM halos 
at the same epoch, from z = to z ~ 10. Our minimal 
approach adopts a physical recipe for star formation with 
dependence on halo assembly time. 

3.1. Halo Assembly Time 

We ad opt the halo a ssemb ly time (redshift) as de- 
fined by lLacev fc Cold (|1993l ) as typical timescale for 
galaxy formation. The assembly redshift z a of a halo 
of mass Mh at redshift z is the redshift at which the 
mass of the main progenitor is which can be 

calc ulated within the extended Prcss-Schcchter formal- 
ism (jBond et al.|[l99lh. For this, we use th e ellipsoidal 
collapse model (|Sheth. Mo fc Tormenl 120011). which re- 
produ ces well numerical simulation results ( Giocoli et al.l 
I2007D . Figure [3] shows typical assembly times at differ- 
ent redshifts in the top-left panel. It is clear that halos 
with the same mass are assembled faster at high z, with 
profound consequences on the UV properties of stellar 
populations, as discussed below. 
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Fig. 2. — Top panel: Star formation rate density (pm* ) an d 
luminosity density (pl) derived by integrating the model UV LFs 
in comparison with the observations before (grey points) and after 
dust correction (black points). Both and pM* are given for 
L > 0.05 L* =3 (M AB < -17.7). Middle: Evolution in the stellar 
mass density (pm+ ) i computed by intergrating the stellar mass 
function to a fixed stellar mass limit of 1O 8 M0. Bottom: Evolution 
in the specific star formation rate (sSFR) as function of rcdshift 
for a galaxy with M, = 5 X 10 9 M Q . 



3.2. Star Formation Modeling 

We populate halos with stars ba sed on the Simple 
Stella r Population (SSP) model s of iBruzual fc Charlotl 
(|2003D . adopting a Salpeter IMF (jSalpeterll 19551) between 
M L = 0.1 M Q and M v = 100 Mq. We use constant 
stellar metallicity Z = 0.02 Z , neglecting redshift evo- 
lution as there is little dependence of the UV luminosity 
on metallicity. We define as l(t) the resulting luminosity 
at 1500 A for a SSP of age t and stellar mass 1 Mq. 

Star formation is assumed to occur in DM halos start- 
ing at their assembly time tH(z a ) with a "burst" contri- 
bution, followed by a constant SFR normalized by l/t age 
("accretion mode"), withi age = tu{z)— tn (z a ) and tn(z) 
is the age of the Universe at redshift z. The physical pic- 
ture is that after an initial burst, each halo accretes the 
same amount of mass and turns it into stars within the 
time span t age . This means that halos of the same mass 
at higher redshifts have higher accretion rates. The re- 
sulting halo luminosity is: 



L(M h , z)=x- [r](M h )M h l{t age )} 



+ (l-x) 



e(M h )M h — / l{t)dt 

*age JQ 



(i) 



The first term r]{Mh) describes the efficiency of the burst 
episode, while the second term e{M h ) describes the ef- 
ficiency of the accretion mode. The free parameter x 
controls the relative contribution of the initial burst to 
the total luminosity of the galaxy at z = 4. 

From Equation [TJ the resulting stellar mass, i.e., the 
time integral of the star formation rate is: 

M*(M h ,z) = x ■ r){M h )M h + (1 - x) ■ e{M h )M h , (2) 

which can be used to obtain stellar densities and specific 
star formation rates. 

3.3. Dust Extinction 

Dust extinction affects significantly the observed UV 
flux, especially at z < 4 (see Figure[2j top panel). Hence, 
it needs to be included i n our mode l . For this, we ap- 
ply dust corrections as in ISmit et al.l(l2012l) : For a spec- 
trum modeled as fx ~ A^, we assume a linear rela- 
tion between the UV-continuum slope j3 and luminosity 

(< ? >= 1iL;( m ^,ab + 19.5) + /3 Muv , AB ; Equation 2 



of lSmit et al.l (j2012f >). Assuming a d ependence of UV ex - 
tinction on /3 as A V v = 4.43+1.99/3 (jMeurer et al.lll999f) . 
and a Gaussian distribution for f3 at each Muv value 
(with dispersion ap — 0.34), then the average < Atjv > 
is given by < A Muv >= 4.43+0 .79 ln(10)crg + 1.99 < (3 > 
(Equation 4 of lSmit et al.ll2012l ). Since negative dust cor- 
rections are non-physical, we adopted the value of for 



any negative < A\jy >. Values for 



d/3 
dM, 



and /3 



are taken from Table 5 of lBouwens et alj ( 2012a|) and are 
listed in the last two columns of our Table[TJ We extrapo- 
lated /?m uv ab to higher and lower redshifts, while letting 
d M vv constant at the z = 4 value, since uncertainities in 
this latter parameter are large. 

3.4. Model Calibration 

To calibrate rj(Mh) and e(M^) we perform abundance 
matching at z = 4, assuming one galaxy per halo and 
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TABLE 1 

Best Fit Schechter Function Parameters for our Model UV LFs and Observed LFs (from the Literature). 





Model Prediction 






Observed LF 




UV continuum parameters 


Rcdshift 


«>*)-3° M* 


a 


{<!>*)-* 


M* a 


Ref 


/3 Muv = -19.5 dMyv 



Z - 


= 0.3 


4.2 ±0.1 


-18.9 ±0.1 


-1.29 ±0.05 


6.2 ± 1.8 


-18.4 ±0.3 


-1.19 ±0.15 


(1) 


-1.45 


-0.13 


Z 


= 1 


1 6+ - 2 
1 - D -o.i 


-19.9 ±0.1 


-1 63+ 04 

1 - oo -0.02 


1.1 ± 0.8 


-20.1 ± 0.5 


-1.63 ±0.45 


(1) 


-1.55 


-0.13 


z 


= 2 


o n+0.2 
z -0.1 


-20.3toi 


-1 60+ 04 
1 - ou -0.06 


2.2 ± 1.8 


-20.2 ± 0.5 


-1.60 ±0.51 


(2) 


-1.70 


-0.13 


z 


= 3 


1.72 ±0.01 


-20.9 + oi 


1 RS+ 005 


1.7 ±0.5 


-21.0 ±0.1 


-1.71 ± 0.53 


(3) 


-1.85 


-0.13 


z 


= 4 


1.30 ±0.01 


-21 0+ ' 2 


i 7 o+0.07 
1 -'°-0.05 


1.3 ±0.2 


-21.0 ±0.1 


-1.73 ±0.05 


(4) 


-2.00 


-0.13 


z 


= 5 


1.4 ±0.1 


-20.6+Q3 


i 77 +0.11 
X - 11 -0.05 


1 4+0-7 
1 '^-0.5 


-20.6 ±0.2 


-1.79 ±0.12 


(5) 


-2.08 


-0.16 


z 


= 6 


1.4 ±0.1 


-20 4+ -' 4 


-1 76+°- 14 
L - ,D -0.12 


1 4+ 11 
-0.6 


-20.4 ±0.3 


-1.73 ±0.20 


(5) 


-2.20 


-0.17 


z 


= 7 


0.9 ±0.1 


-20.2 ±0.2 


-1 84+° 12 


QQ + 0.7 


-20.1 ± 0.3 


-2.01 ± 0.21 


(6) 


-2.27 


-0.21 


z 


= 8 


0.5 ±0.1 


-20 2+ 4 
zu - z -0.2 


-1 Q 2 + 011 
^ -0.15 


4+ 04 


— 20 s+ ' 3 


-l.98±g;l 


(7) 


-2.34 


-0.25 


z ■ 


= 10 


0.2 ±0.1 


-19.74±g;g 


9 1O+0.25 
z - lo -0.02 


0.1 ± 0.1 


-19.6 6 


-1.73 6 


(8) 


C 


C 



References . (1) lArnouts et"aL1 J2005I"); (2) IQesch et all l20Tol). (3) IReddv fc Steidell (t2009Tl ; (4) IBouwens et~aTI J2007D ; (5) 

IBouwens et~aTl [|2012bl) ; f 6) IBouwens et al.1 l|201H) ; C7) IBradlev et al.1 1 <2012f l; (81 IQesch et all 1120121) . 

Note. — The best fit Schechter Function parameters, as a function of redshift (column 1), for our model-predicted LFs (columns 2-4; 
see Figure [TJ, and for observed UV LFs (taken from the literature; columns 5-8). Quoted errors for our model predictions are derived by 
propagating the uncertainty in the L(M^,2; = 4) calibration (see also Figure(3]l. The last two columns show the adopted slopes (3m uv =- 19.5 

and intercepts — to the UV-continuum slope /3 to UV luminosity relationship as in IBouwens etaLI <|2012al . Table 5). 
a Units: 10- 3 Mpc" 3 . 

b Values were fixed for the fit of the Schechter function. 
c No dust correction at z ~ 10. 



equating the number of galaxies with luminosity greater 
than L (after dust correction) to the number of halos 
with mass great than M^: 

/ n(M h ,z = 4)dM h = / (f>(L,z = 4)dL, (3) 

where n(M^ z) is the mass func t ion o f DM halos ob- 
tained adopting ISheth fc Tormenl (|1999t ) mass function. 
This gives us a luminosity versus halo mass relation at 
z — 4, L(Mh, z — 4), shown in the bottom-right panel of 
Figure |3] From this we can then infer rj(Mh) and e(Mh) 
by solving Equation [T] (bottom-left panel of Figure |3]). 
We calibrate these two quantities independently, so that 
their linear combination also satisfies L(Mh, z = 4) by 
construction. The shaded areas in both panels represent 
the uncertainty in the model calibration derived by vary- 
ing the z — 4 LF par ameters within the la confidence 
regions in Figure 3 of IBouwens et al.l (|2007| ). From the 
bottom- left panel of Figure [3] it is immediate to see that 
halos with Mh ~ 10 11 — 10 12 M Q have the highest specific 
star formation efficiencies. This is not surprising, given 
the shapes of the LFs and DM mass functions, and it may 
reflect negative feedback from Supernovae at the low- 
mass end, and dissipationless merging and AGN feed- 
back at the high-mass end (e.g.. lCoorav fc Milosavljevic I 
[2005h . 

Our final calibration step is selecting a value for the 
only free parameter in the model, x, i.e., the contribution 
of the burst to the total luminosity at z = 4. To this 
purpose, we compute model predictions over the redshift 
range < z < 8 with varying x values, and adopt the 
value of x which minimizes the residuals relative to the 
observed LFs. The best match to observations is given 
by x = 0.1, that is a 10% of contribution from the initial 
burst to the total halo luminosity at z = 4; the model is 
however not very sensitive to the exact value for as long 



as x -€i 1 (see top panels of Figure 01 ■ 

4. RESULTS AND DISCUSSION 

The bottom-right panel of Figure contains our pre- 
dictions for the observed L(Mh, z) relation at different 
redshifts. A decreasing contribution from dust is the 
main cause for the brightening of the relation towards 
higher redshifts at fixed halo mass. The predictions of 
our model for the LFs over the z < 10 time span are 
shown overplotted to the observations in Figure [TJ The 
model reproduces both low-z and high-z UV LFs remark- 
ably well, suggesting that the evolution of the UV LF 
across most of cosmic time can indeed result from the 
redshift evolution of the DM halo mass function, cou- 
pled with simple star formation histories beginning at 
the halo assembly time. 

The model, calibrated to the observed Schechter LF at 
z = 4, produces Schechter functions at all other epochs. 
Furthermore, the predicted LFs well approximates the 
Schechter functions with the observed best-fit parame- 
ters reported in Table Q] (Figure [lj. In particular, the 
model correctly describes the evolution of the faint-end 
slope a, from its shallow low-z value a(z ~ 0) ~ —1.3 to 
the steepening observed at z > 6, where a < — 1.7. At 
high-z, the model L Fs are similar to the LFs predicted by 
iTrenti et al.1 ((2010) ; in addition, our new model is also 
successful in reproducing the observed LFs also at low 
redshifts, all the way down to z ~ 0, resolv ing the puz- 
zling quick rise of a from z ~ to 1 (e.g., IQesch et al.l 
12011 . 

Figure [2] shows the model predictions for the red- 
shift evolution of the star formation rate density pM t (z) 
and luminosity density Pl(z) (top panel). The model- 
observations agreement for pm„ (z) is again very good at 
all epochs from z w to z w 8. At z ~ 1 the model 
appea rs to over-predict pm* as measured by IQesch et al.l 
( 2012) by ~ 0.8 dex. This might be due to very short as- 
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Fig. 3. — Upper-left panel: DM halo assembly time t age versus halo mass and redshift comput ed on the basis of the Extended Press- 
Schechter formalism. Upper-right panel: Calibration of our model using the Bouwcns et al. (2007) LF at z = 4. Bottom-left panel: Star 
formation efficiencies ^(M^) and e(Mj,) for the burst and accretion mode, respectively. Halo masses in the range M h ~ lO^-lO 12 Mq have 
the highest efficiencies. Bottom-right panel: The relation between galaxy luminosity and DM halo mass, L(Mf l ,z), plotted at z = 2,4,8. 
The main driver of LF brightening with increasing redshift at fixed halo mass is the decrease in dust content and younger ages. In all 
panels, shaded areas represent 68% confidence regions; these are derived by propagating the uncertainty in the L(M^, z = 4) calibration. 



sembly times for z ~ 10 halos (t age < 10 8 Myr), hence to 
a dominant contribution to the UV light from the very 
young stellar populations produced in the burst mode. 
Another possibility is sample variance in the observa- 
tions: The LF at thi s frontier redshift is based on a 
single LBG detection (jBouwens et al.lf2011b lOesch et all 
|2012|) . and may need revision. In fact. iZheng et al.ll2012l 
derive from a gravitationally lensed source in CLASH 

PmA z = !0) = (1.8lf;?) x 1CT 3 M Q Mpc- 3 yr _ \ in full 
agreement with our model predictions. 

Figure |4] (Panels (a)-(b)) shows model-observations 
comparisons for the z = 2 and 6 LFs; the model results 
are plotted for different values of the burst fractional con- 
tribution x to the total UV luminosity, and for differ- 
ent dust extinction corrections (as para metrized by /?). 
Varying /3 within the uncertainty given bv lBouwens et all 
(|2012af) has little impact on the predicted LFs, which 
are thus robust against uncertainties in the amount and 



treatment of dust obscuration. In contrast, the predicted 
LFs do depend on the choice of x. As the star forma- 
tion histories are modeled to approach a single episode 
of star formation at the halo assembly time, i.e., x — > 1, 
the L(Mh, z) relation shows an increasingly stronger de- 
pendence on the assembly time: at low redshift, halos 
become too faint relative to the observations, and the 
LFs and pl are under-estimated; at high redshifts, halos 
are too UV-bright and LFs and pl are overestimated. 

An outcome of the model is a very modest contribution 
of the initial burst phase to the UV luminosity of galax- 
ies at all epochs. The left-bottom axes of panel (c) of 
Figure U] show the contribution of the burst mode to the 
total galaxy luminosity, for halos of different masses. The 
burst contribution to the galaxy luminosity becomes in- 
creasingly less important towards lower redshifts, as t age 
increases. The redshift evolution of the L{, urst / L to t ratio 
is faster for smaller halos, as for these t age evolves faster 
(Figure [3]). The initial star formation burst contributes 
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Fig. 4. — (a)-(b): Impact on the 2 = 2 and 2 = 6 LFs given by varying the relative contribution of burst to galaxy luminosity. Values 
x <C 1 give overall a good description of the observed LFs, while x ~ 1 underestimates the LF at low z. Additionally in panel (a), the 
impact of variation of /3m uv =— 19 5 by ±30% is shown as grey area, (c): Contribution of the burst mode to the total luminosity as a 
function of redshift and halo mass (left-bottom axis) and to the total stellar mass as a function of halo mass, since there is no redshift 
dependence (right-top axis), (d): Stellar masses as a function of UV luminosity for 2 ~ 4 and 7 (NC: nebular contamination corrected data 
with equivalent width evolution (EW(z)) and without (EW), respectively.) 



however at least as much as the slow-accretion phase 
to the stellar mass. Figure 21 Panel (d), shows stellar 
mass as a function of magnitude M^ jr.nn- The re d 
points show the relation at z — 4 of iStark et al.l {2009), 
and the blue points are re-normalized for accounting for 
emission lines (see IStark et al.l (|2012f )). We are slightly 
below the z — 4 observations, but at z — 7 we pre- 
dict M+ — M^B.is nn relation which matches the data of 
IStark et al.l <|20ia i. corrected for emission line contami- 
nation and with an equivalent width of nebular emission 
increasing with redshift. 

At fixed halo mass, Equation[5]does not depend on red- 
shift. The contribution of the burst phase to the stellar 
mass is roughly 50% at Mh ~ 10 14 Af Q , and increases to 
about 70% at Mh ~ 1O 9 M (see up-right axes of panel 
(c) of Figure |4j . While the UV bright stars produced 



in the burst rapidly fade away, the burst, happening by 
construction at the halo assembly time, quickly sets the 
stellar mass of the luminous galaxy associated to a given 
halo to roughly half of the final mass. This is consistent 
with a picture in which star formation rate in luminous 
galaxies closely follows the DM halo accretion rate (e.g., 
iPeng et al.|[2010HFaucher-Giguere et al.ll201lD . 

The middle Panel of Figure [2] shows the comparison 
between model and observations for pm, ■ Overall, the 
model fits the data well at all redshifts. At z = 7, we 
under-estimate by 0.2 5 dex the (not emission-corrected) 
iGonzalez et al.l (|2011h data point. On the other hand , 
we are broadly consistent with the IStark et al.l ([2012D 
measurements. The sSSFR is shown in the bottom panel 
of Figure [2] and shows overall a good agreement with the 
observations. 
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5. CONCLUSION 

We have presented a model for the evolution of the 
UV LF based on the simple assumption that all mas- 
sive DM halos host a galaxy with a star formation his- 
tory that is closely related to the halo assembly time. 
We specifically adopt for the star formation histories the 
combination of an initial star formation burst, a proxy 
for a short-phase, intense star formation event, and a 
constant sSFR, a proxy for secular, low-level star forma- 
tion activity fueled by steady gas accretion. While the 
assumption that each DM halo hosts only one galaxy 



is clearly a simplification, especially towards lower red- 
shifts, the model is remarkably successful in reproducing 
major features of the evolving star forming galaxy pop- 
ulation since z ~ 8; this is also thanks to the fact that 
massive galaxies which, at later times, will share a com- 
mon halo, will be mostly quenched of their star formation 
activity. It is remarkable that our simple parametriza- 
tion reproduces very well the evolution of the UV LFs 
over the whole cosmic time since z ~ 8 down to z ~ 0, as 
well as the evolution of the cosmic specific star formation 
rate, luminosity density and stellar mass density. This 
clearly demonstrates the key role played by DM halos in 
shaping the properties of their galaxies. 
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